Edge-dependent selection rules in magic triangular graphene flakes 
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The electronic shell and supershell structure of triangular graphene quantum dots has been studied 
using density functional and tight-binding methods. The density functional calculations demonstrate 
that the electronic structure close to the Fermi energy is correctly described with a simple tight- 
binding model where only the Pz orbitals perpendicular to the graphene layer are included. The 
results show that (i) both at the bottom and at the top of the pz band a supershell structure similar 
to that of free electrons confined in a triangular cavity is seen, (ii) close to the Fermi level the shell 
structure is that of free massless particles, (iii) triangles with armchair edges show an additional 
sequence of levels ('ghost states') absent for triangles with zigzag edges while the latter exhibit edge 
states, and (iv) the observed shell structure is rather insensitive to the edge roughness. 
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Recent experimental success in manufacturing single 
layer graphene flakes on various surfaces [TJ HJ [31 2] has 
made graphene a new playground for theoretical and 
computational physics, [3 [6l [7l 13 [9] and more and more 
experimental results are emerging. [TOl [TTJ [12] Most of 
the recent interest has been focused in the effects caused 
by the peculiar band structure of graphite near the Fermi 
level (ef): Electrons and holes behave as massless parti- 
cles (Dirac fermions) due to the linear dispersion relation 
although their velocity is very small. [13] 

The triangular shape of two-dimensional clusters is 
particularly interesting because, in the case of free elec- 
trons, it supports perhaps the most persistent and regular 
supershell structure of all systems. '14] Furthermore, the 
triangular shape is preferred in two-dimensional metallic 
systems [Tallin] and in plasma clusters. [T7] For tetrava- 
lent elements, triangular clusters have been observed in 
silicon. [18j It is reasonable to expect that such shapes 
can be observed also for carbon, and this is supported 
further by the fact that equilateral triangles of graphene 
can be cut with the two most stable edge structures, the 
zigzag edge and the armchair edge. 

In this letter, we wish to point out that finite graphene 
flakes (or quantum dots) have an intriguing energy spec- 
trum close to the Fermi level. We have performed 
electronic structure calculations for triangular graphene 
flakes using the density functional theory (DFT) for all 
the valence electrons, and a tight-binding (TB) approach 
that considers only the carbon pz electrons (Hiickel 
model) . Our results show that already in small triangular 
flakes (A^=300, L=5 nm) the electronic levels close to ep 
can be understood as those of free massless electrons con- 
fined in a triangular cavity. Especially, we demonstrate 
that the edge structure has a selective role in the elec- 
tronic shell structure: The zigzag edge prohibits a whole 
sequence of localized states inside the cluster although 
it supports edge states. This leads to well-defined edge- 
dependent selection rules that are based on an analytical 
model. Recently, Yamamoto et al. addressed the pres- 



ence (absence) of edge states at ep in zigzag (armchair) 
triangles of graphene, and the effect on the optical ab- 
sorption, |19j but the simple principles of the underlying 
energy spectrum have remained unexplained. 
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FIG. 1: Cross-over of the valence and conduction bands at 
the Fermi energy (top) , and the density of states (bottom, pz 
electrons) of an infinite graphene sheet. 

It is well-known that the atomic Pz electrons perpen- 
dicular to the graphene plane are responsible for the cap- 
tivating band structure shown in Fig. [l] with the va- 
lence the conduction bands meeting at the corners of the 
hexagonal Brillouin zone. [501111] The Fermi surface con- 
sists of a discrete set of these points of high-A: value, and 
the resulting density of states (DOS) has a zero weight 
at ep- The cross-over regions have locally hourglass-like 
shapes, which results in the linear and isotropic electron 
dispersion relation in the conduction band (e > ei?=0) 
but only in a small energy interval. Since the atomic p^ 
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electrons are perpendicular to the graphene plane their 
interaction with the neighboring atoms does not have 
any directional dependence and, consequently, they can 
be described as s-type electrons in the TB model. Ne- 
glecting also the differential overlap between atomic sites 
the system can be described with the traditional Hiickel 
model 



—t, if i,j nearest 
0, otherwise 



neighbours 



where the hopping parameter t (resonance integral) de- 
termines the width of the bands and the on-site energy 
is chosen to be ep = 0. We choose to present our results 
in units t = 1. The resulting TB bands (Fig. [T| reaches 
from -3 to +3 (in real graphene our unit t corresponds to 
about 2.6 eV). 

A conceptual cutting of a finite graphene flake breaks 
covalent bonds yielding edges with dangling bonds. We 
consider the dangling bonds to be passivated, say with 
hydrogen. Since the covalent bonding with hydrogen in- 
volves sp^ hybridized orbitals, the passivation is expected 
to have only a small effect on the perpendicular pz elec- 
tron states. Therefore, we neglect this effect in our TB 
model and follow Areshkin et al. [2T] and treat the edge 
atoms in the same footing as bulk atoms. Moreover, we 
will completely neglect the interaction of graphene with 
the possible substrate and treat the graphene flake as an 
isolated two-dimensional cluster or quantum dot. As we 
shall see, the results of the simple TB model agree well 
with those of the full DFT calculation. 
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FIG. 2: Upper panel: DOS at the bottom of the TB band 
shown as a function oi q = y/e + 3t. Blue: zigzag triangle with 
10000 atoms; red: armchair triangle with 9918 atoms; and 
black: result for free electrons in a triangular cavity. Lower 
panel: Electron densities of the six lowest energy levels. 

It has been shown that at the bottom of the valence 



band the TB model gives exactly the free electron states 
for a triangular lattice, [22] and the same is true also for 
the hexagonal graphene. Consequently, at the bottom 
(and at the top) the energy levels are expected to show 
the same shell structure as free electrons in a triangular 
cavity, determined by the equation [531 121] 



eo(n^ 



(1) 



where ep = Sn^h^ /3meL^ , L being the length of the tri- 
angle side. The quantum numbers must satisfy m > 1 
and n > 2m. Determination of the electron effective 
mass in the graphene lattice gives for the TB model 
eo = 47r^</9iV, where N is the number of atoms in the tri- 
angle {L — 3d\/N/2 for a large triangle, d is the nearest- 
neighbor distance). 

The shell structure manifests itself as a regular varia- 
tion of the density of the states (DOS) which can be de- 
termined by Gaussian convolution of the discrete levels. 
Figure [2] shows DOS close to the bottom of the valence 
band obtained from the above equation and compared to 
the TB model for two graphene triangles, one with 10000 
atoms (zigzag edge) and the other with 9918 atoms (arm- 
chair edge). The proflles are clearly similar and exhibit 
the beating pattern of the supershcU structure [17j . Note 
that DOS is plotted as a function of + 3t making the 
shells equidistant. Fig. [2] shows also the electron den- 
sities corresponding to the six lowest energy levels (for 
degenerate states we show the sum of the density). The 
density patterns are identical to those of free electrons 
conflned in a triangle [16) or wave modes in triangular 
resonators 

The Fermi level of graphene consists of two equivalent 
points at the border of the Brillouin zone (see Fig. [ij 
where the conduction and valence bands open as circular 
cones resulting a linear dispersion relation for electrons 
e(k) — Chk, where C is the velocity. Thus, it is to be 
expected that the electron dynamics is not determined 
by the Schrodinger equation but by the wave equation 
of massless particles (or the Dirac equation). For free 
particles conflned in a triangle the energy levels are still 
determined by Eq. ([T|) but now it results in the square 
of the energy, i.e 



n'^ + m'^ — nm, 



(2) 



where ei = 2TTt/^/3N. It is interesting to note that these 
energy levels were actually computed for the wave equa- 
tion much earlier than for the Schrodinger equation. |23j 
Figure [3] shows TB-DOS above the Fermi energy for 
two large triangles (~ 10000 atoms) with zigzag and arm- 
chair edges and compares them with the levels of free 
massless electrons [Eq. ([2|]. The results are the follow- 
ing: (i) Each energy level has an additional degeneracy 
of two due to the two equivalent points at e^. (ii) The 
zigzag triangle shows the levels of Eq. Q with index 
values m > 1 and n > 2m while the armchair edge shows 
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FIG. 3: Upper panel: TB-DOS at the Fermi level displayed 
as a function of energy (red and blue curves) compared to the 
density of levels of Eq. ([2| (black curves). The zigzag triangle 
has 10000 atoms and the armchair triangle 9918 atoms. Lower 
panel: DOS of the full DFT calculation for the triangular 
C321H51 (zigzag, red) and CsaoHeo (armchair, blue) flakes. 
The inset shows the levels above the Fermi surface where the 
zigzag spectrum is shaded. 



all the levels where n > m > 1. (iii) The states are 
much less dense than at the bottom of the band, and 
Eq. ^ describes only the lowest states accurately, (iv) 
Due to the sparseness of the states, no supershell oscilla- 
tions are visible for the massless particles (although the 
supershell structure of ordinary electrons is clearly seen 
m Fig. [2]). (v) The zigzag edge supports particularly vis- 
ible edge states {27, ,28 that appear at ei? as a prominent 
peak. The number of these states equals the number of 
the outermost edge atoms in zigzag triangles, which is 



In order to compare our results with a more realistic 
calculation, we have performed DFT calculations for tri- 
angular C321H51 (zigzag) and C330H60 (armchair) flakes 
with the CPMD program. [23] The DFT calculations use 
a plane wave basis set {Ecut = 50 Ry), pseudopotentials, 
[30] and a generalized gradient-corrected PBE approx- 
imation for the exchange-correlation energy. [31] The 
resulting DFT-DOS of all valence electrons is plotted in 
Fig. |3jb) for both systems, and they show overall features 
characteristic for graphite. The zigzag edge states at 
are visible, and the closest conduction states obey the 
simple analytical model of Eq. ([2]). The even- numbered 
peaks are split for the armchair triangle, which is a result 
reproduced by TB (the splitting reduces with increasing 
system size). 



FIG. 4: Electron density of the 3rd [(a) and (b)] and 5th 
[(c) and (d)] energy levels above the Fermi energy in arm- 
chair triangles ('ghost states', labeled in Fig. [3| each has a 
degeneracy two), (a) and (c) are computed for a large TB 
triangle of 4920 C atoms while (b) and (d) are DFT results 
for a C330H60 molecule. 



The lowest conduction states that are labeled in Fig. 
[3] show fascinating details, and the electron densities of 
two such states are visualized in Figure [4| For compar- 
ison, we show the same states/orbitals calculated for a 
large triangle with the TB model (4920 C atoms) and for 
a small triangle calculated with the DFT method (330 C 
atoms) . The internal structure (symmetry) of the states 
is clearly similar, and therefore, it is independent of the 
triangle size and the model used. The states close to the 
Fermi level appear very different from those at the bot- 
tom of the band (Fig. |2|. They are no i simple densities of 
massless particles confined in a triangle since the density 
profile does not decay to zero at the edges. The corre- 
sponding electron levels are close to the Brillouin zone 
boundary, having large fc-values, and the wave functions 
have pronounced oscillations with wave lengths that are 
related to the unit cell size. These oscillations guarantee 
that the wave function will be formally zero at the edges, 
but the corresponding pseudowave function of the mass- 
less particle does not necessarily show the same behavior. 
An interesting feature in Fig. [4] is that the states have 
simple geometric structure of triangular symmetry. The 
size (number) of the triangles decreases (increases) with 
increasing energy, i.e. the pattern repeats itself. These 
'ghost states' are completely absent for the zigzag trian- 
gles, and they correspond to quantum numbers of Eq. 
([2]) not allowed for free electrons in a triangular box [i.e. 
2m > n > m > 1 in Eq. (|2l]. 
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FIG. 5: Electron density (TB model) of the 2nd [(a) and (b)] 
and 4th [(c) and (d)] energy levels above the Fermi energy 
(labeled in Fig. [3| for armchair and zigzag triangles of 4920 
and 5181 C atoms, respectively. 



Figure |5] shows the electron densities corresponding to 
the 'normal' low energy states that obey the standard 
selection rules (to > 1 and n > 2m). Again, the electron 
density does not necessarily vanish at the edges of the 
triangle. The corresponding states for the armchair and 
zigzag triangles display obvious differences despite the 
fact that they involve the same set of quantum numbers 
(and energy). 

Finally, we want to note that a small roughness of the 
edge does not remove the peculiar states shown in Fig. 
|4]or change the shell structure close to the Fermi level. 
These 'ghost states' form a triangular network, and it 
would be interesting to study if they can exist also in the 
graphene flakes with hexagonal, parallelogram, or trape- 
zoidal shapes. 

In conclusion, we have computed the electronic struc- 
ture of triangular graphene flakes and shown that the 
DOS proflle close to ep is independent of the triangle 
size, and it can be described with the simple TB model. 
The zigzag flakes exhibit well-known edge states, and the 
armchair triangles show an additional set of 'ghost states' 
(different selection rules) where the corresponding elec- 
tron density makes a triangular lattice. In large triangles 
of 5000-10000 C atoms, the energy levels can be described 
accurately by considering free massless particles confined 
in a triangular cavity. Presumably, the electronic states 
near the Fermi surface are not sensitive to the dielec- 
tric substrate, and we expect that these fascinating wave 
functions can be observed with scanning tunneling mi- 
croscopy. 
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